IMPACT  OF  CHANNEL  STATISTICS  AND  CORRELATION  ON  UNDERWATER 
ACOUSTIC  COMMUNICATION  SYSTEMS 

by 

JESSE  SCOTT  CROSS 

A  THESIS 

Presented  to  the  Faculty  of  the  Graduate  School  of  the 
MISSOURI  UNIVERSITY  OF  SCIENCE  AND  TECHNOLOGY 
in  Partial  Fulfillment  of  the  Requirements  for  the  Degree 

MASTER  OF  SCIENCE  IN  ELECTRICAL  ENGINEERING 

2011 

Approved  by 


Dr.  Yahong  Rosa  Zheng 
Dr.  Chengshan  Xiao 
Dr.  Steven  Grant 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

2011 


2.  REPORT  TYPE 


4.  TITLE  AND  SUBTITLE 

Impact  of  Channel  Statistics  and  Correlation  on  Underwater  Acoustic 
Communication  Systems 

6.  AUTHOR(S) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Missouri  University  of  Science  and  Technology, 1870  Miner 
Circle, Rolla, MO, 65409 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


3.  DATES  COVERED 

00-00-2011  to  00-00-2011 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

Several  statistical  properties  of  underwater  acoustic  channels  gathered  from  experiment  data  are  analyzed. 
The  baseband  channel  impulse  response  (CIR)  is  estimated  using  a  time  domain  least  squares  technique 
with  a  sliding  window  applied  to  the  probing  sequences.  From  the  CIR  estimation,  the  probability 
distribution  functions  (PDFs)  of  the  magnitude,  real  part,  imaginary  part,  and  phase  of  the  CIR  are 
calculated.  Gamma  Rayleigh,  and  compound  k  distributions  are  fitted  to  the  magnitude  PDF  and  the 
fitness  of  the  distributions  are  calculated  with  a  two-sample  Kolmogorov-Smirnov  test.  Other  statistics 
such  as  the  autocorrelation  function,  coherence  time,  and  scattering  function  are  evaluated.  The  results 
show  that  the  underwater  acoustics  channels  are  worse  than  the  Rayleigh  fading  commonly  seen  as  the 
worst  case  radio  channel.  Furthermore,  the  spatial  and  intertap  correlation  matrices  of  multiple  input 
multiple  output  (MIMO)  systems  are  estimated  using  experimental  data.  It  is  shown  that  underwater 
acoustic  MIMO  channels  exhibit  high  spatial  and  temporal  correlation.  The  bit  error  rate  (BER)  of  the 
receiver  using  Frequency-domain  turbo  equalization  is  also  evaluated  in  different  channel  correlation 
setups,  demonstrating  strong  effects  of  the  spatial-temporal  correlation  function  on  the  performance. 

15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF 

ABSTRACT 

OF  PAGES 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

58 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


Ill 


PUBLICATION  THESIS  OPTION 


This  thesis  contains  papers  published  in  the  style  used  by  the  Navy  Journal  of  Un¬ 
derwater  Acoustics  (NJUA)  and  in  the  style  used  by  the  IEEE  Military  Communications 
Conference  (MILCOM).  Pages  3-19  have  been  published  in  the  NJUA  in  2010.  Pages  20- 
32  have  been  published  in  MILCOM  2010  and  pages  33-45  are  accepted  for  publication 
in  MILCOM  2011. 


ABSTRACT 


Several  statistical  properties  of  underwater  acoustic  channels  gathered  from  exper¬ 
iment  data  are  analyzed.  The  baseband  channel  impulse  response  (CIR)  is  estimated 
using  a  time  domain  least  squares  technique  with  a  sliding  window  applied  to  the  prob¬ 
ing  sequences.  From  the  CIR  estimation,  the  probability  distribution  functions  (PDFs)  of 
the  magnitude,  real  part,  imaginary  part,  and  phase  of  the  CIR  are  calculated.  Gamma, 
Rayleigh,  and  compound  k  distributions  are  fitted  to  the  magnitude  PDF  and  the  fitness 
of  the  distributions  are  calculated  with  a  two-sample  Kolmogorov- Smirnov  test.  Other 
statistics  such  as  the  autocorrelation  function,  coherence  time,  and  scattering  function 
are  evaluated.  The  results  show  that  the  underwater  acoustics  channels  are  worse  than 
the  Rayleigh  fading  commonly  seen  as  the  worst  case  radio  channel. 

Furthermore,  the  spatial  and  intertap  correlation  matrices  of  multiple  input  mul¬ 
tiple  output  (MIMO)  systems  are  estimated  using  experimental  data.  It  is  shown  that 
underwater  acoustic  MIMO  channels  exhibit  high  spatial  and  temporal  correlation.  The 
bit  error  rate  (BER)  of  the  receiver  using  Frequency- domain  turbo  equalization  is  also 
evaluated  in  different  channel  correlation  setups,  demonstrating  strong  effects  of  the 
spatial-temporal  correlation  function  on  the  performance. 


V 


ACKNOWLEDGMENTS 


I  would  like  to  extend  my  gratitude  toward  all  who  have  helped  me  achieve  my 
master’s  degree. 

First,  I  would  like  to  thank  my  mother  and  all  the  rest  of  my  family  members  for 
their  support  while  I  worked  toward  my  master’s  degree. 

I  would  like  to  thank  Dr.  Yahong  Rosa  Zheng  for  her  guidance  throughout  my  time 
at  MS&T,  for  the  undergraduate  student  research  position  from  June  2007  to  May  2010, 
and  for  the  graduate  research  position  starting  from  June  2010  to  Dec.  2011.  I  would 
also  like  to  thank  Dr.  Chengshan  Xiao  and  Dr.  Steven  Grant  for  being  a  part  of  my 
master’s  degree  committee.  1  would  like  to  thank  Jian  Zhang  for  his  help  with  the  papers 
included  in  this  thesis. 

1  would  like  to  express  my  gratitude  for  the  chancellor’s  fellowship.  I  am  grateful 
for  the  support  provided  by  the  Intelligent  Systems  Center.  The  papers  included  in  this 
thesis  were  supported  in  part  by  the  National  Science  Foundation  under  Grant  ECCS- 
0846486  and  the  Office  of  Naval  Research  under  Grants  N00014-10- 1-0174. 


VI 


TABLE  OF  CONTENTS 


Page 

PUBLICATION  THESIS  OPTION .  iii 

ABSTRACT  .  iv 

ACKNOWLEDGMENTS  .  v 

LIST  OF  ILLUSTRATIONS .  ix 

LIST  OF  TABLES .  x 

SECTION 

1.  INTRODUCTION .  1 

PAPER 

1.  Statistical  Properties  of  Multiple-Input  Multiple-Output  (MIMO)  Underwater 

Acoustic  Communications  Channels .  3 

Abstract .  3 

1.1.  Introduction .  3 

1.2.  Channel  model  and  channel  statistics .  5 

1.2.1.  Channel  Estimation .  6 

1.2.2.  Magnitude  Distributions .  7 

1.2.3.  Kolmogorov- Smirnov  (KS)  Test .  8 

1.2.4.  Channel  Scattering  Function  .  8 

1.3.  Experimental  results .  9 

1.4.  Conclusion .  17 

1.5.  Acknowledgment  .  18 

References .  18 

2.  Statistical  Channel  Modeling  of  Wireless  Shallow  Water  Acoustic 

Communications  from  Experiment  Data .  20 


vii 

Abstract .  20 

2.1.  Introduction .  20 

2.2.  Channel  model  and  channel  statistics .  21 

2.2.1.  Channel  Estimation .  22 

2.2.2.  Magnitude  Distributions .  23 

2.2.3.  Kolmogorov- Smirnov  (KS)  Test .  24 

2.2.4.  Channel  Scattering  Function  .  24 

2.3.  Experimental  results .  25 

2.4.  Conclusion .  30 

2.5.  Acknowledgment  .  31 

References .  31 

3.  Impact  of  Spatial  Correlation  of  Fading  Channels  on  the  Performance  of  MIMO 

Underwater  Communication  Systems .  33 

Abstract .  33 

3.1.  Introduction .  33 

3.2.  MIMO  Acoustic  Communication  Systems .  34 

3.2.1.  Channel  Estimation .  35 

3.2.2.  Frequency  Domain  Turbo  Equalization .  35 

3.2.3.  Spatial  Correlation .  37 

3.3.  Experiment  Results .  38 

3.3.1.  Spatial  Correlation .  39 

3.3.2.  Transceiver  BER  Performance .  41 

3.4.  Simulation  Results  .  42 

3.5.  Conclusion .  44 

3.6.  Acknowledgment  .  44 

References .  44 


Vlll 


SECTION 

2.  CONCLUSIONS .  46 

VITA  .  47 


IX 


LIST  OF  ILLUSTRATIONS 

Figure  Page 

1.1  Estimated  time- varying  frequency-selective  channel  impulse  responses  in  the 

RACE08  experiment .  10 

1.2  Estimated  channel  impulse  responses  in  ACOMM09  experiment .  11 

1.3  The  PDFs  of  the  RACE08  complex  CIR .  12 

1.4  The  PDFs  of  the  ACOMM09  complex  CIR .  13 

1.5  The  normalized  autocorrelation  of  the  real  part,  imaginary  part,  and  complex 

envelop  of  the  CIR .  15 

1.6  The  normalized  cross  correlation  between  the  real  and  imaginary  parts  of 

the  CIR .  16 

1.7  The  normalized  spatial  correlation  between  Rxl  and  the  other  hydrophones 

for  RACE08  channel .  16 

1.8  The  channel  scattering  function  of  the  RACE08  experiment .  17 

2.1  Estimated  time- varying  frequency- selective  channel  impulse  responses.  .  .  27 

2.2  The  PDFs  of  the  complex  CIR .  28 

2.3  The  normalized  autocorrelation  of  the  real  part,  imaginary  part,  and  complex 

envelop  of  the  CIR .  29 

2.4  The  normalized  cross  correlation  between  the  real  and  imaginary  parts  of 

the  CIR .  30 

2.5  The  channel  scattering  function .  30 

3.1  Magnitude  of  the  channel  correlation  matrix  of  a  2x2  MIMO  channel  ....  39 

3.2  Spatial  correlations  between  different  receivers .  40 

3.3  Simulation  results .  43 


X 


LIST  OF  TABLES 

Table  Page 

1.1  The  first  and  second  moments  of  Rayleigh,  Gamma,  and  Compound  KDis- 

tributions .  8 

1.2  The  Kolmogorov-Smirnov  results  for  the  measured  PDF  when  compared  to 

different  fitted  distributions .  14 

2.1  The  first  and  second  moments  of  Rayleigh,  Gamma,  and  Compound  k  Dis¬ 
tributions  .  24 

2.2  The  Kolmogorov-Smirnov  results  for  the  measured  PDF  when  compared  to 

different  fitted  distributions .  28 

3.1  4 >Tx  and  for  hydrophones  at  200  m .  41 

3.2  Transmission  distance  and  receiver  spacing  BER  results .  42 


1.  INTRODUCTION 


Underwater  acoustic  (UWA)  channels  have  a  much  longer  range  (1  km  -  1000  km) 
than  underwater  radio  channels;  however,  underwater  acoustic  channels  are  much  more 
severe  and  present  many  additional  challenges  such  as  excessive  path  loss,  a  complicated 
propagation  paths,  and  a  time-varying  nature.  Statistical  analysis  of  in-air  radio  fre¬ 
quency  (RF)  channels  is  already  established.  In  the  past,  most  modeling  of  underwater 
acoustic  channels  has  focused  on  wave  propagation  with  limited  statistical  analysis.  Re¬ 
cent  ocean  experiments  have  shown  that  the  UWA  channels  exhibit  worse  than  Rayleigh 
fading  and  an  extreme  Doppler  shift  in  comparison  to  the  carrier  frequency.  The  long  mul¬ 
tipath  delay  time  results  in  severe  intersymbol  interference  (ISI)  and  a  channel  impulse 
response  (CIR)  that  has  many  more  taps  than  an  RF  system.  All  of  these  properties  re¬ 
sult  in  the  need  to  re-estimate  the  channel  repeatedly  which  increases  the  computational 
complexity  of  an  underwater  communication  system  as  compared  to  an  RF  system. 

The  papers  presented  in  this  thesis  analyze  data  gathered  in  three  ocean  experi¬ 
ments.  The  first  is  Reschedule  Acoustic  Communication  Experiment  (RACE)  conducted 
in  Narragansett  Bay,  Rhode  Island,  March  2008;  the  second  is  the  ACOMM  experiment 
conducted  in  New  Jersey  in  2009  (ACOMM09);  and  the  third  is  the  Surface  Process  and 
Acoustic  Communications  Experiment  conducted  in  fall  2008  (SPACE08)  by  the  Woods 
Hole  Oceanographic  Institute  (WHOI).  Data  from  these  experiments  are  used  to  inves¬ 
tigate  the  statistical  characteristics  of  UWA  channels,  including  the  probability  density 
function  (PDF),  autocorrelation  function,  crosscorrelation  function,  and  scattering  func¬ 
tion  of  the  baseband  equivalent  CIRs.  A  sliding  window  least  squares  technique  is  used 
to  estimate  and  track  the  time- varying  CIRs.  Various  theoretical  PDFs  are  fitted  to  the 
estimated  CIR  PDF  and  the  Kolmogorov-Smirnov  test  is  used  to  quantize  the  level  of 
fitness.  The  effects  of  spatial  correlation  on  the  bit  error  rate  (BER)  are  investigated. 
Trials  with  differing  amounts  of  receivers  but  with  equal  aperture  are  run.  A  simula¬ 
tion  is  conducted  to  confirm  the  experimental  results  are  due  to  spatial  correlation.  The 
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Kronecker  properties  of  the  spatial  correlation  matrix  are  analyzed.  The  matrix  is  de¬ 
composed  into  a  series  of  Kronecker  products  and  the  error  resultant  matrix  formed  by 
the  decomposed  matrices  and  the  original  matrices  is  quantified. 

ft  is  shown  that  the  experimental  PDF  matches  the  compound  k  distribution.  The 
compound  k  distribution  presents  additional  challenges  over  the  Rayleigh  distribution 
of  in  air  RF  channels.  Since  the  taps  follow  the  compound  k  distribution  rather  than 
the  Rayleigh  distribution,  there  is  a  much  larger  probability  that  the  channel  tap  will 
be  low  while  at  the  same  time  there  is  an  equal  probability  that  the  tap  will  be  large. 
This  necessitates  a  receiver  with  a  much  larger  dynamic  response  so  that  the  receiver 
has  a  good  response  with  low  tap  values  but  does  not  become  saturated  with  high  tap 
values.  The  channel  coherence  time  of  the  channels  analyzed  is  shown  to  be  about  470 
symbol  periods  long.  This  implies  that  the  receiver  must  use  data  blocks  smaller  than 
this  in  order  to  have  an  invariant  channel  during  the  transmission  of  the  data  block. 
When  the  number  of  receivers  is  reduced  while  maintaining  a  constant  aperture,  the 
BER  remains  relatively  constant.  This  shows  spatial  correlation  plays  a  large  part  in 
receiver  placement. 
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PAPER 

1.  Statistical  Properties  of  Multiple-Input  Multiple-Output  (MIMO)  Underwater 

Acoustic  Communications  Channels 
Jesse  Cross1,  Jian  Zhang2,  and  Yaliong  Rosa  Zheng1 
department  of  Electrical  and  Computer  Engineering, 

Missouri  University  of  Science  and  Technology,  Rolla,  MO  65409 
2  now  at  Broadcom  Corp.,  San  Diego,  CA  92127 


Abstract 

Statistical  characteristics  of  underwater  acoustic  channels  are  estimated  by  ocean 
experimental  data  obtained  in  the  RACE08  and  ACOMM09  experiments.  First,  base¬ 
band  complex  channel  impulse  responses  (CIRs)  are  estimated  by  a  time-domain  least 
square  method  with  sliding  windows  applied  to  long  probing  sequences.  Second,  the 
statistics  of  the  CIRs  are  analyzed,  including  the  probability  density  functions  (PDF)  of 
the  real  and  imaginary  parts  and  of  the  magnitude  and  phase,  autocorrelation  functions, 
channel  coherence  time,  and  scattering  function.  Spatial  correlation  of  multiple  trans¬ 
ducers  and  hydrophones  are  also  extracted  for  multiple-input  multiple-output  (MIMO) 
channels. 

1.1.  Introduction 

Underwater  acoustic  (UWA)  propagation  is  proven  to  be  the  effective  means  of 
underwater  wireless  communication  for  medium  and  long  ranges  (1  km  -  1000  km). 
However,  underwater  acoustic  channels,  especially  shallow  water  horizontal  channels,  are 
often  more  challenging  than  radio  frequency  (RF)  channels  due  to  excessive  path  loss, 
complicated  propagation  environment,  and  time- varying  inhomogeneous  media.  Statis¬ 
tical  modeling  of  RF  channels  has  been  well  established  in  the  literature  and  has  played 
important  roles  in  capacity  analysis,  transceiver  design,  and  performance  evaluation. 


4 


However,  analysis  and  modeling  of  UWA  channels  focused  on  the  acoustic  wave  propa¬ 
gation  mechanism  in  the  past  and  limited  studies  are  available  for  extracting  statistical 
properties  [1,2]. 

Extensive  ocean  experiments  have  been  conducted  over  the  past  two  decades  and 
existing  UWA  channel  statistics  studies  show  that  UWA  channels  are  often  worse  than 
Rayleigh  fading,  which  in  turn  heavily  impacts  the  bit  error  rate  (BER)  performance  [4]. 
Multipath  and  Doppler  effect  of  UWA  are  described  in  [1,  2]  providing  guidelines  in 
equalizer  design.  Wideband  UWA  channel  is  modeled  in  [5]  recently.  It  is  shown  that 
UWA  channel  capacity  is  also  greatly  impacted  by  accurate  channel  models  [6]. 

This  paper  analyzes  ocean  experiment  data  collected  in  the  Reschedule  Acous¬ 
tic  Communication  Experiment  (RACE),  Narragansett  Bay,  Rhode  Island,  March  2008, 
and  the  data  collected  in  the  ACOMM09  experiment  conducted  in  New  Jersey  in  2009 
(ACOMM09).  This  paper  also  investigates  the  statistical  characteristics  of  UWA  chan¬ 
nels,  including  the  probability  density  function  (PDF),  autocorrelation  function,  cross- 
correlation  function,  and  scattering  function  of  the  baseband  equivalent  channel  impulse 
responses  (CIRs).  The  complex  CIRs  are  estimated  by  a  sliding  window  least  square 
method,  where  small  window  length  and  sliding  width  are  used  for  long  probing  data 
sequences  to  accurately  track  the  time-varying  CIRs.  The  Kolmogorov- Smirnov  (KS) 
test  [7]  is  employed  to  statistically  test  the  fitness  between  the  PDF  of  estimated  CIRs 
and  theoretical  PDFs.  The  test  results  show  that  the  experimental  channel  magnitude 
PDF  matches  the  compound  K  distribution  with  small  shape  parameters.  This  results 
complement  previous  investigation  in  wireless  RF  [8, 9]  and  underwater  [4, 10]  channels. 
Our  results  also  show  that  the  phases  of  CIRs  are  approximately  uniform  in  [0,  27r]  and 
the  autocorrelation  function  exhibits  exponential  decay  that  is  dominated  by  the  Gamma 
component  of  the  compound  K  distribution.  The  channel  scattering  function  is  obtained 
from  the  correlation  function  providing  channel  characteristics  in  the  delay-Doppler  plane. 
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These  statistical  properties  provide  design  guidelines  for  channel  estimation,  channel 
equalization,  pilot  insertion,  and  transceiver  signaling. 

1.2.  Channel  model  and  channel  statistics 

Consider  first  a  single-input  single-output  communication  system  with  one  trans¬ 
ducer  and  one  hydrophone.  Let  s(t)  denote  the  baseband  equivalent  complex  signal  and 
S(t)  =  Re{s(i)eJ'2,r^ct}  denote  the  passband  transmit  signal  with  a  carrier  frequency  fc. 
At  the  receiver,  the  received  signal  is  R(t)  =  R e{f(£)eJ'27I'ct},  where  r(t)  is  baseband 
complex  envelop.  Then  r(t)  can  be  represented  as 

r(t)  —  J  h(t,r)s(t  —  r)dr  (1.1) 

where  h(t,  r)  is  the  equivalent  complex  baseband  impulse  response  for  a  time-varying 
frequency-selective  channel.  In  a  multipath  and  Doppler  spread  propagation  environ¬ 
ment,  we  assume  there  are  /  multi-paths  and  each  path  has  distinguished  delay  and 
Doppler  shift.  Then,  the  channel  impulse  response  is  modeled  as 

i 

hit ,  r)  =  Ai  ■  ( -,2 r;  fen)  .  ^  ^  (1.2) 

i=  1 

where  At,  rir  and  fda  are  the  gain,  propagation  delay,  and  instantaneous  Doppler  shift  at 
the  i-th  path.  Although  acoustic  communication  channels  are  time  varying  resulting  in 
time- varying  Aj,  t*,  and  fd.% ,  the  time  variation  is  very  small  within  the  channel  coherence 
time,  thus  these  parameters  are  usually  assumed  to  be  constant  in  a  short  time  interval. 
Therefore  we  drop  the  dependence  of  the  channel  impulse  response  with  respect  to  t  and 
denote  h{r)  as  the  so  called  quasi-static  channel. 

Next,  consider  a  MIMO  communication  system  with  Nt  transducers  and  Nr  hy¬ 
drophones.  Denote  hm,n{l)  the  subchannel  CIR  linking  the  m-th  hydrophone  with  the 
n-th  transducer  with  sampling  at  the  symbol  duration  Ts ,  where  /  =  1,  2,  •  •  •  ,  L.  The 
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MIMO  channel  is  represented  by  the  channel  matrix 

hi,i(l)  h\p{l)  •••  h\ ,Nt(l) 

h(o=  :  :•••:,  (1.3) 

hi, i(l)  hi, 2(1)  ■■■  h]\jr,Nt  (0 

and  the  input-output  relationship  of  the  MIMO  system  is  then 

L 

r(k)  —  y H(l)s(k  —  l)  +  w(fc),  (1.4) 

1=1 

where  k  is  the  time  index,  transmitted,  received,  and  noise  vectors  are  respectively  s(k)  = 
[si(k)  ■■■sNt(k))T,  r (k)  =  [f^k)  ■■■rNr(k)]T,  and  w (k)  =  [w^k)  ■■■wNr(k)]T,  with 
superscript  T  denoting  transpose. 

1.2.1.  Channel  Estimation.  A  time-domain  least  squares  (LS)  method  is 
employed  in  this  paper  to  probe  the  channel  impulse  response  (CIR)  from  the  baseband 
input  and  output  signals.  A  long  probing  sequence  is  transmitted  by  the  transducer  and 
the  CIRs  are  estimated  progressively  by  using  a  sliding  window  of  size  Np. 

Let  denote  the  probing  symbols  of  the  n-tli  transducer  at  the  p-th  step.  The 
estimated  time-domain  CIR  for  the  m-th  hydrophone  is  =  [(h)),  , )T,  •  •  •  ,  (h^  Ni)T]T, 
where  hm,n  =  [hsmn{  1),  •  •  •  ,  hsrnn(L)]T  with  a  channel  length  L  and  Nt  is  the  number  of 
transducers,  can  be  estimated  as 

K  =  K)'  •  C 


(1.5) 
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where  ()*  denotes  the  pseudo  inverse,  ,r^n{Np)]T  which  is  the  corre- 


sponding  MIMO  channel  output  vector,  and  sp  =  | 

Tfl-'-ITS,,]  with 

<*(£) 

1) 

<(1) 

TAP  — 
n 

<(L  +  1) 

<(L)  ■  ■ 

<(2) 

.  <{NP) 

<(Nr-  1) 

'  7  +  1)_ 

where  n  —  1,  •  •  •  ,  Nt.  Let  Nd  denote  the  sliding  step  of  the  data  window.  The  CIR  is 
re-estimated  every  NdTs  interval,  where  Ts  is  the  symbol  period. 

1.2.2.  Magnitude  Distributions.  Let  X  have  a  Rayleigh  distribution  with  a 
parameter  a ,  then  the  conditional  PDF  of  X  is  given  as 

/*(*)  =  p(-^l)  (1-7) 

The  PDF  of  a  Gamma  distributed  A"  is  given  as 

^  =  <l8> 

where  r(-)  is  the  Euler  gamma  function  and  a  and  b  are  scalar  parameters. 

The  unconditioned  PDF  of  X,  if  X  has  a  compound  K  distribution,  is  given  as 


fx(x) 


4  ( x  yK 


(1.9) 


where  Kv-i  is  the  modified  bessel  function  of  the  second  kind  of  order  v  —  1.  A  Rayleigh 
distribution  results  when  v  — >  oo  and  (3v  =  2 a'2  remains  constant.  Therefore,  as  v 
decreases  the  channel  gets  worse  than  Rayleigh  fading. 

The  first  and  second  moments  were  used  to  fit  the  PDFs  of  the  different  distributions 
to  the  measured  data  and  are  given  in  the  Table  1.1. 


Table  1.1  The  first  and  second  moments  of  Rayleigh,  Gamma,  and  Compound  KDistri- 
butions  _ 


Rayleigh 

Gamma 

Compound  K 

E(x) 

__  Hk 
aV  2 

ab 

~ 737  r(i/+o.5) 

2  r(^) 

E(x 2) 

I 

2  U 

ab2 

(3v 

1.2.3.  Kolmogorov-Smirnov  (KS)  Test.  The  KS  test  is  one  of  the  most  useful 
and  general  nonparametric  tools  to  quantify  the  similarity  of  two  empirical  distribution 
functions.  The  two-sample  KS  test  can  serve  as  a  goodness  of  fit  test,  and  is  used  to 
determine  if  two  datasets  share  a  common  probability  distribution.  No  assumptions  are 
made  about  the  distributions  of  the  two  datasets.  The  null  hypothesis  of  the  KS  test  is 
that  both  datasets  do  not  differ  significantly  and  therefore  follow  the  same  distribution. 
The  alternative  hypothesis  is  that  the  two  datasets  follow  different  distributions.  The 
Kolmogorov-Smirnov  statistic  D  is  defined  as  [7] 


Dny  =  sup  | Fn(x)  -  Fn>(x) | 

X 


(1.10) 


where  Fn(x)  and  Fn>(x )  are  the  empirical  cumulative  distribution  functions  (CDF)  of  the 
two  distributions  to  be  compared  and  sup  is  the  supremum  of  the  set  \Fn(x)  —  Fn>(x)  \. 
The  alternative  hypothesis  will  be  accepted  at  level  a  if 


Qa  ^  F)r 


nn 


n  +  n' 


(1.11) 


where  Qa  is  found  by  P(Q  <  Qa )  =  1  —  a  and  Q  is  a  Kolmogorov  distributed  random 
variable;  n  and  n'  are  the  number  of  samples  in  the  compared  distributions,  respectively. 

1.2.4.  Channel  Scattering  Function.  Generally,  the  channel  impulse  response 
h(t,  t )  is  considered  as  a  two-dimensional  random  process  for  frequency  selective  channels. 
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The  general  second-order  autocorrelation  function  of  h(t,  r)  can  be  written  as 

Rh{tl,t2',T1,T2)  =  E[h(t1,T1)h*(t2,T2)].  (1.12) 

where  E()  denotes  statistical  expectation.  To  simplify  this  four-dimensional  correlation 
function,  we  make  an  assumption  that  the  channel  is  wide  sense  stationary  in  the  short 
term,  i.  e,  the  second-order  statistics  depends  only  on  the  difference  between  time  instants, 
rather  than  the  absolute  time.  Therefore,  the  autocorrelation  of  h(t,  r)  is  simplified  as 

E (h(t  +  At,  ti)  •  h*(t,  r2))  =  Rh(At,  ti)<5(ti  -  r2)  (1.13) 

and  the  channel  scattering  function  is  defined  by  [11] 

/OO 

Rh(At,r)e-j2*fAtdAt  (1.14) 

-oo 

The  two-dimensional  scattering  function  (1.14)  depicts  the  power  spectral  density  of 
the  channels  along  geo-time  and  completely  describes  the  second-order  statistics  of  a 
stationary  fading  channel. 

1.3.  Experimental  results 

The  underwater  experiment,  named  as  Reschedule  Acoustic  Communication  Ex¬ 
periment  (RACE),  was  conducted  in  Narragansett  Bay,  Rhode  Island,  by  Woods  Hole 
Oceanographic  Institution  (WHOI),  in  March  2008.  The  water  depth  varied  between  9 
to  14  meters.  The  transmitter  with  an  array  of  two  transducers  was  fixed  on  a  tripod  and 
was  four  meters  above  the  sea  bottom.  The  receiver  with  an  array  of  12  hydrophones  was 
also  fixed  on  a  tripod  and  was  located  two  meters  above  the  sea  bottom.  In  this  study, 
the  data  was  collected  when  the  distance  between  the  transmitter  and  receiver  was  1000 
meters.  The  carrier  frequency  was  fc  =  11.5  kHz,  the  sampling  rate  fs  =  39.0625  kHz, 
and  the  bandwidth  B  =  fsj 10  =  3.90625  kHz. 
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A  probing  sequence  with  a  duration  of  22  seconds  was  transmitted  to  probe  the 
channel  impulse  responses.  A  sliding  window  was  used  with  a  window  size  Np  =  128  and 
the  interval  between  adjacent  window  =  32.  Hence,  the  sampling  rate  at  the  absolute 
time  domain  is  1/(32  •  Ts)  =  122  Hz  which  means  the  channels  are  estimated  in  every  8.2 
milliseconds.  The  channel  impulse  responses  (CIRs)  were  estimated  progressively  by  the 
LS  method  described  in  (1.5).  For  each  tap,  2144  samples  were  obtained  in  one  probing 
data  packet.  The  two-dimensional  time-varying  CIRs  for  the  subchannel  corresponding 
to  (Txl,  Rxl)  and  (Tx2,  Rxl)  are  plotted  in  Fig.  1.1(a)  and  Fig.  1.1(b).  The  frequency- 
selective  channel  length  spans  25  symbol  periods. 


Txl - Rxl  Tx2 - Rxl 


(a)  Estimated  CIRs  for  the  subchannel  Txl — (b)  Estimated  CIRs  for  the  subchannel  Tx2 — 
Rxl  Rxl 


Figure  1.1  Estimated  time- varying  frequency-selective  channel  impulse  responses  in  the 
RACE08  experiment. 


The  ACOMM09  was  conducted  by  Naval  Research  Laboratory  (NRL)  at  the  coast¬ 
line  of  New  Jersey  in  May  2009.  The  symbol  period  was  0.2  ms  and  the  carrier  frequency 
was  17  kHz.  Modulations  included  QPSK  (Quadrature  Phase  Shift  Keying),  8-PSK, 
and  16-QAM  (Quadrature  Amplitude  Modulation).  The  transmit  array  consisted  of  two 
transducers  located  at  25  ~  75  m  deep  from  the  sea  surface.  Two  receive  hydrophone 
arrays,  ACDS2  and  ACDS3,  were  deployed,  with  their  distances  to  the  transmitter  being 
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2  km  and  3  km,  respectively.  The  depths  of  ACDS2  and  ACDS3  were  20  ~  34.4  m  and 
40  ~  54.4  m,  respectively.  Each  hydrophone  array  consisted  of  eight  hydrophones  with 
the  inter- hydrophone  spacing  being  2.06  m. 

The  channel  length  was  measured  as  L  =  200  in  this  experiment.  An  example  of  the 
estimated  channels  for  two-transducer  transmission,  as  shown  in  Fig.  1.2,  demonstrated 
that  the  two  transmitter  channels  were  unbalanced  with  Tx2  having  approximately  half 
of  the  power  as  those  of  Txl.  The  CIR  of  Tx2-Rx4,  for  example,  was  also  non-minimum 
phase  which  may  cause  significant  difficulties  in  channel  synchronization  and  equalization. 


Txl — Rx3  Tx2 — Rx3 


Txl — Rx5  Tx2 — Rx4 


Figure  1.2  Estimated  channel  impulse  responses  in  ACOMM09  experiment. 


Figures  1.3(a)  and  1.3(b)  show  the  real  and  imaginary  parts  of  the  RACE08  mea¬ 
sured  data  distribution  versus  a  Gaussian  distribution  with  the  same  mean  and  different 


12 


variances.  Fig.  1.4(a)  and  Fig.  1.4(b)  show  the  measurements  for  ACOMM09.  We  can 
see  the  measured  PDFs  of  the  real  and  imaginary  parts  do  not  follow  the  Gaussian  dis¬ 
tribution.  For  the  probability  densities  around  large  values,  the  Gaussian  distribution 
with  large  variance  matches  the  tails  of  the  measured  PDFs.  However,  for  the  probability 
densities  at  small  values,  the  Gaussian  distributions  with  small  variance  fit  the  peaks  of 
the  measured  PDFs  better. 


Normalized  real  part  Normalized  imaginary  part 


(a)  Real  Part  PDF  (b)  Imaginary  Part  PDF 

Magnitude  PDF  Phase  PDF 


(c)  Magnitude  PDF 


(d)  Phase  PDF 


Figure  1.3  The  PDFs  of  the  RACE08  complex  CIR.  (a)  PDF  of  real  part  versus  Gaussian 
with  a  =  1.6  and  0.75.  (b)  PDF  of  imaginary  part  versus  Gaussian  with 
a  =  1.25  and  0.66.  (c)  PDF  of  magnitude  versus  compound  Kand  other 
distributions,  (d)  PDF  of  phase  versus  uniform  distribution. 
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Normalized  Real  Part  Normalized  Imaginary  Part 


(a)  Real  Part  PDF 


(b)  Imaginary  Part  PDF 


Magnitude  PDF 


Normalized  Amplitude 


Phase  PDF 


(c)  Magnitude  PDF 


(d)  Phase  PDF 


Figure  1.4  The  PDFs  of  the  ACOMM09  complex  CIR.  (a)  PDF  of  real  part  versus 
Gaussian  with  a  =  0.55  and  0.40.  (b)  PDF  of  imaginary  part  versus  Gaussian 
with  a  =  0.50  and  0.39.  (c)  PDF  of  magnitude  versus  compound  Kand  other 
distributions,  (d)  PDF  of  phase  versus  uniform  distribution. 


Since  the  real  and  imaginary  parts  are  not  well  represented  by  Gaussian  distribu¬ 
tions,  it  is  reasonable  to  predict  the  magnitude  probability  distribution  will  not  follow  a 
Rayleigh  distribution.  Fig.  1.3(c)  shows  the  RACE08  measured  magnitude  PDF  versus 
fitted  compound  Kand  other  distributions.  Fig.  1.4(c)  shows  the  same  for  the  ACOMM09 
experiment.  It  can  be  seen  visually  that  the  Rayleigh  and  Gamma  distributions  are  not 
a  good  fit.  The  KS  test  was  used  to  compare  the  measured  PDF  with  different  distri¬ 
butions  to  find  the  best  match  to  the  measured  PDF.  The  KS  test  results  for  different 
distributions  are  depicted  in  Table  1.2.  The  p-value  represents  the  probability  of  a  test 
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Table  1.2  The  Kolmogorov-Smirnov  results  for  the  measured  PDF  when  compared  to 
_ different  fitted  distributions _ 


Gamma 

Rayleigh 

Compound  K 

RACE08 

Parameters 

p-value 

Dn,n’ 

Results 

a  =  1.065,  b  =  0.8 
1.18  x  10"6 

0.16 

reject 

a  =  0.70 
1.1  x  10"51 
0.48 
reject 

u  =  0.63,  (3  —  6.04 
0.96 

0.07 

accept 

ACOMM09 

Parameters 

p-value 

Dn,n 

Results 

a  =  1.05,  b  =  0.4 
1.10  x  10"7 
0.7646 
reject 

a  =  0.4 

1.0  x  10"50 
0.7546 
reject 

u  =  1.25,  /3  —  0.42 
0.1196 

0.0392 

accept 

statistic  at  extreme  as  the  one  actually  observed.  If  the  significance  level  a  =  5%,  the 
null  hypothesis  is  rejected  when  the  p-value  is  less  than  0.05.  The  Dn  is  the  maximum 
distance  between  two  CDFs,  and  the  results  have  two  choices:  reject  or  accept  the  null 
hypothesis.  The  fit  test  parameters  for  different  distributions  were  chosen  from  the  KS 
test  that  returned  the  highest  p-value.  For  the  RACE08  data,  the  null  hypothesis  that 
the  measured  distribution  follows  the  fitted  compound  Kdistribution  with  v  =  0.63  and 
/3  =  6.04  was  accepted  with  a  p-value  of  0.96  and  Dn  n>  =  0.07.  However,  the  null  hypoth¬ 
esis  that  the  measured  data  is  subject  to  the  best  fit  Rayleigh  distribution  was  rejected 
with  a  very  small  p-value,  1.1  x  10~51.  The  null  hypothesis  that  the  measured  channel 
followed  the  fitted  Gamma  distribution  was  rejected  with  a  p-value  of  1.18  x  10-6.  There¬ 
fore,  we  conclude  that  the  measured  data  follows  a  compound  Kdistribution  rather  than 
a  Rayleigh  or  Gamma  distribution  for  the  R.ACE08  channel.  Similarly,  the  ACOMM09 
measured  channel  also  followed  a  compound  K  distribution  with  v  =  1.25  and  (3  =  0.42, 
as  shown  in  Table  1.2. 

The  normalized  autocorrelation  functions  of  each  channel  tap  for  real  part,  imag¬ 
inary  part  and  complex  envelop  are  depicted  in  Fig.  1.5.  It  is  noted  that  these  three 
autocorrelation  curves  behave  similarly  all  with  exponential  decay.  Define  the  channel 
coherence  time  as  the  time  over  which  the  autocorrelation  is  above  0.5.  The  coherence 
time  was  approximately  Tc  =  0.12s  for  this  channel.  It  implies  that  the  channel  can 
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Figure  1.5  The  normalized  autocorrelation  of  the  real  part,  imaginary  part,  and  com¬ 
plex  envelop  of  the  CIR.  The  autocorrelation  function  was  dominated  by  the 
Gamma  component  of  the  CIR.  The  coherence  time  of  the  channel  was  about 
120  ms. 


be  considered  to  be  invariant  if  the  length  of  data  block  is  less  than  Tc  ■  B  ps  470  sym¬ 
bols.  The  relation  between  the  coherence  time  and  the  maximum  Doppler  spread  is  given 


by  [12] 
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16  rcfm 


(1.15) 


where  fm  is  maximum  Doppler  shift.  In  this  case,  fm  =  1.5  Hz  and  the  dynamic  range  of 
the  instantaneous  Doppler  spread  is  [—1.5, 1.5]  Hz.  Figure  1.6  shows  the  cross  correlation 
of  the  real  and  imaginary  part  of  the  CIR.  It  is  seen  that  the  real  parts  and  the  imaginary 
part  are  uncorrelated. 

The  spatial  correlation  between  hydrophone  Rxl  and  the  other  hydrophones  is 
presented  in  Fig.  1.7  for  the  RACE08  channel,  where  Rx  1  was  the  hydrophone  located 
at  the  very  bottom  of  the  vertical  array  and  Rxl2  was  the  one  at  the  top  of  the  array. 
Define  the  coherence  distance  as  the  distance  away  from  Rxl  where  the  spatial  correlation 
is  above  0.5.  It  is  clear  from  the  figure  that  within  the  channel  coherence  time  120ms,  the 
spatial  correlation  between  Rxl  and  Rx2  was  above  0.5,  thus  the  two  receive  hydrophones 
were  spaced  within  the  coherence  distance.  The  spatial  correlation  between  Rx  1  and 
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Crosscorrelation  of  real  and  imaginary  parts 


Figure  1.6  The  normalized  cross  correlation  between  the  real  and  imaginary  parts  of  the 
CIR.  The  real  and  imaginary  parts  were  uncorrelated. 


Figure  1.7  The  normalized  spatial  correlation  between  Rxl  and  the  other  hydrophones 
for  RACE08  channel.  The  normalized  correlation  was  averaged  across  the  25 
taps. 


other  hydrophones  were  all  below  0.5  and  can  be  considered  as  spatially  uncorrelated. 
The  spatial  correlation  between  other  pairs  of  hydrophones  exhibited  similar  results  as 
being  correlated  if  spaced  within  coherent  distance  and  uncorrelated  if  spaced  far  apart. 
The  figures  are  not  shown  here  for  brevity. 


17 


The  two-dimensional  scattering  function  is  shown  in  Fig.  1.8  by  performing  Fourier 
transform  on  the  correlation  function.  We  see  that  the  average  Doppler  shift  is  very  small, 
close  to  zero  Hz  since  the  transceiver  platforms  were  always  fixed  during  the  experiment 
and  the  motion  of  the  wave  was  not  significant.  From  the  power  spectral  density,  the 
channel  has  a  Doppler  spread  of  [—1.5, 1.5]  Hz  at  a  delay  of  1  ~  2.5  ms. 


Scattering  function 


-8  -6  -4  -2  0  2  4  6  8 

Doppler  Spread  (Hz) 


Figure  1.8  The  channel  scattering  function  of  the  RACE08  experiment. 

1.4.  Conclusion 

In  this  paper,  we  presented  the  statistical  characteristics,  including  PDF  and  second- 
order  statistics,  for  the  UWA  channels  measured  in  RACE08  and  ACOMM09  ocean  ex¬ 
periments.  Experimental  results  demonstrate  that  the  real  and  imaginary  parts  of  the 
complex  UWA  channels  do  not  follow  Gaussian;  the  PDF  of  the  magnitude  fits  the  com¬ 
pound  K  distribution;  the  PDF  of  phase  approximately  follows  uniform  distribution. 
The  autocorrelations  of  CIRs  behave  like  decaying  exponential  curves  and  the  channel 
coherence  time  is  also  roughly  estimated.  The  very  low  spatial  correlation  between  the 
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receivers  has  a  profound  impact  on  the  design  of  MIMO  systems.  The  scattering  function 
provides  better  understanding  on  the  spectral  density  of  the  UWA  channels. 
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Abstract 

This  paper  analyzes  statistical  characteristics  of  underwater  acoustic  channels  es¬ 
timated  by  ocean  experimental  data.  The  baseband  complex  channel  impulse  responses 
(CIRs)  are  estimated  by  a  time-domain  least  square  method  with  sliding  windows  ap¬ 
plied  to  long  probing  sequences.  The  probability  density  functions  (PDF)  of  the  real 
part,  imaginary  part,  magnitude,  and  phase  of  the  CIR  are  evaluated,  and  the  two- 
sample  Kolmogorov-Smirnov  test  is  used  to  measure  fitness  of  the  magnitude  PDF  to 
the  Gamma,  Rayleigh,  and  compound-K  distributions.  The  second-order  statistics  of 
the  channel  are  also  investigated  in  terms  of  autocorrelation  function,  channel  coherence 
time,  and  scattering  function.  The  experimental  results  demonstrate  that  underwater 
channels  are  often  worse  than  Rayleigh  fading  channels. 

2.1.  Introduction 

Underwater  acoustic  (UWA)  propagation  is  proven  to  be  the  effective  means  of 
underwater  wireless  communication  for  medium  and  long  ranges  (1  km  -  1000  km). 
However,  underwater  acoustic  channels,  especially  shallow  water  horizontal  channels,  are 
often  more  challenging  than  radio  frequency  (RF)  channels  due  to  excessive  path  loss, 
complicated  propagation  mechanism,  and  time-varying  inhomogeneous  media.  Statisti¬ 
cal  modeling  of  RF  channels  has  been  well  established  in  the  literature  and  has  played 
important  roles  in  capacity  analysis,  transceiver  design,  and  performance  evaluation. 
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However,  analysis  and  modeling  of  UWA  channels  focused  on  the  acoustic  wave  propa¬ 
gation  mechanism  in  the  past  and  limited  studies  are  available  for  extracting  statistical 
properties  [1,2]. 

Extensive  ocean  experiments  have  been  conducted  over  the  past  two  decades  and 
existing  UWA  channel  statistics  studies  show  that  UWA  channels  are  often  worse  than 
Rayleigh  fading,  which  in  turn  heavily  impacts  the  bit  error  rate  (BER)  performance  [4,4], 
Multipath  and  Doppler  effect  of  UWA  are  described  in  [1,  2]  providing  guidelines  in 
equalizer  design.  Wideband  UWA  channel  is  modeled  in  [5]  recently.  It  is  shown  that 
UWA  channel  capacity  is  also  greatly  impacted  by  accurate  channel  models  [6]. 

This  paper  analyzes  ocean  experiment  data  collected  in  the  Reschedule  Acoustic 
Communication  Experiment  (RACE),  Narragansett  Bay,  Rhode  Island,  March  2008,  and 
investigates  the  statistical  characteristics  of  UWA  channels,  including  the  probability  den¬ 
sity  function  (PDF),  autocorrelation  function,  crosscorrelation  function,  and  scattering 
function  of  the  baseband  equivalent  channel  impulse  responses  (CIRs).  The  complex 
CIRs  are  estimated  by  a  sliding  window  least  square  method,  where  small  window  length 
and  sliding  width  are  used  for  long  probing  data  sequences  to  accurately  track  the  time- 
varying  CIRs.  The  Kohnogorov-Smirnov  (KS)  test  [7]  is  employed  to  statistically  test 
the  fitness  between  the  PDF  of  estimated  CIRs  and  theoretical  PDFs.  The  test  results 
show  that  the  experimental  channel  magnitude  PDF  matches  the  compound  K  distri¬ 
bution  with  small  shape  parameters.  Our  results  also  show  that  the  phases  of  CIRs  are 
approximately  uniform  in  [0,  2i r]  and  the  autocorrelation  function  exhibits  exponential 
decay  that  is  dominated  by  the  Gamma  component  of  the  compound-K  distribution.  The 
channel  scattering  function  is  obtained  from  the  correlation  function  providing  channel 
characteristics  in  the  delay-doppler  plane. 

2.2.  Channel  model  and  channel  statistics 

Let  s(t)  denote  the  baseband  equivalent  complex  signal  and  S(t)  =  R e{s(t)eJ'27r^ct} 
denote  the  passband  transmit  signal  with  a  carrier  frequency  fc.  At  the  receiver,  the 
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received  signal  is  R(t)  =  Re{r(t)e-?27r-'ct},  where  r(t)  is  baseband  complex  envelop.  Then 
r(t)  can  be  represented  as 

r(t)  —  J  h(t,r)s(t  —  r)dr  (2.1) 

where  h(t,  r)  is  the  equivalent  complex  baseband  impulse  response  for  a  time-varying 
frequency-selective  channel.  In  a  multipath  and  Doppler  spread  propagation  environ¬ 
ment,  we  assume  there  are  /  multi-paths  and  each  path  has  distinguished  delay  and 
Doppler  shift.  Then,  the  channel  impulse  response  is  modeled  as 

i 

h(t,  t)  =  A  ■  (yi::{LlJ  fcTi)  ■  5(t  -  n )  (2.2) 

i—  1 

where  Ai,  r*,  and  fd,%  are  the  gain,  propagation  delay,  and  instantaneous  Doppler  shift  at 
the  i-th  path.  Although  the  time- varying  channel  results  in  time- varying  Ai,  r,;,  and  fd.n 
these  parameters  can  be  assumed  to  be  constant  in  a  short  time  interval,  which  is  the  so 
called  quasi-static  channel. 

2.2.1.  Channel  Estimation.  A  time-domain  least  squares  (LS)  method  is 
employed  in  this  paper  to  probe  the  channel  impulse  response  (CIR)  from  the  baseband 
input  and  output  signals.  A  long  probing  sequence  is  transmitted  by  the  transducer  and 
the  CIRs  are  estimated  progressively  by  using  a  sliding  window  of  size  Np. 

Let  si  denote  the  probing  symbols  of  the  n-th  transducer  at  the  p- th  step.  The 
estimated  time-domain  CIR  for  the  m-th  hydrophone  is  hfn  =  [(hfn  l)T,  •  •  •  ,  (h^  ArJr]T, 
where  hm>n  =  [h^  ( 1),  •  •  •  ,  hsmn(L)]T  with  a  channel  length  L  and  Nt  is  the  number  of 
transducers,  can  be  estimated  as 
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where  ()*  denotes  the  pseudo  inverse,  ,r^n{Np)]T  which  is  the  corre¬ 

sponding  MIMO  channel  output  vector,  and  sp  =  [Tf|  •  •  •  \TpNt]  with 
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<(L~  1)  '■ 

<(1) 

TP— 

<(£+!) 

<(L)  ■  ■ 

*S(  2) 
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'  sn(^p“^+l) 

(2.4) 


where  n  —  1,  •  •  •  ,  Nt.  Let  Nd  denote  the  sliding  step  of  the  data  window.  The  CIR  is 
re-estimated  every  NdTs  interval,  where  Ts  is  the  symbol  period. 

2.2.2.  Magnitude  Distributions.  Let  A"  have  a  Rayleigh  distribution  with  a 
parameter  a ,  then  the  conditional  PDF  of  X  is  given  as 


fx(x) 


X 


( 7 “ 


exp( 


(2.5) 


The  PDF  of  a  Gamma  distributed  X  is  given  as 

<2-6> 

where  T(-)  is  the  Euler  gamma  function  and  a  and  b  are  scalar  parameters. 

The  unconditioned  PDF  of  X,  if  X  has  a  compound  K  distribution,  is  given  as 


fx(x) 


4  ( x  yK 


(2.7) 


where  Kv_\  is  the  modified  bessel  function  of  the  second  kind  of  order  v  —  1.  A  Rayleigh 
distribution  results  when  v  — >  oo  and  f3v  =  2<x2  remains  constant.  Therefore,  as  v 
decreases  the  channel  gets  worse  than  Rayleigh  fading. 

The  first  and  second  moments  were  used  to  fit  the  PDFs  of  the  different  distributions 
to  the  measured  data  and  are  given  in  the  Table  2.1. 
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Table  2.1  The  first  and  second  moments  of  Rayleigh.  Gamma,  and  Compound  k  Distri¬ 
butions  _ 


Rayleigh 

Gamma 

Compound  K 

E(x) 

ab 

7/^  r(i/+o.5) 

2  r(^) 

E(x 2) 

W  ' 

ab2 

(3v 

2.2.3.  Kolmogorov-Smirnov  (KS)  Test.  The  KS  test  is  one  of  the  most  useful 
and  general  nonparametric  tools  to  quantify  the  similarity  of  two  empirical  distribution 
functions.  The  two-sample  KS  test  can  serve  as  a  goodness  of  fit  test,  and  is  used  to 
determine  if  two  datasets  share  a  common  probability  distribution.  No  assumptions  are 
made  about  the  distributions  of  the  two  datasets.  The  null  hypothesis  of  the  KS  test  is 
that  both  datasets  do  not  differ  significantly  and  therefore  follow  the  same  distribution. 
The  alternative  hypothesis  is  that  the  two  datasets  follow  different  distributions.  The 
Kolmogorov-Smirnov  statistic  D  is  defined  as 


Dny  =  sup  | Fn(x)  -  Fn>(x) | 

X 


(2.8) 


where  Fn(x)  and  Fn>(x )  are  the  empirical  cumulative  distribution  functions  (CDF)  of  the 
two  distributions  to  be  compared  and  sup  is  the  supremum  of  the  set  \Fn(x)  —  Fn>(x)  \. 
The  alternative  hypothesis  will  be  accepted  at  level  a  if 


Qa  ^  Dr 


nn 


n  +  n' 


(2.9) 


where  Qa  is  found  by  P(Q  <  Qa)  =  1  —  a  and  Q  is  a  Kolmogorov  distributed  random 
variable;  n  and  n'  are  the  number  of  samples  in  the  compared  distributions,  respectively. 

2.2.4.  Channel  Scattering  Function.  Generally,  the  channel  impulse  response 
h(t,  t )  is  considered  as  a  two-dimensional  random  process  for  frequency  selective  channels. 
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The  general  second-order  autocorrelation  function  of  h(t,  r)  can  be  written  as 

Rh{ti,t2',T1,T2)  =  E[h(t1,T1)h*(t2,T2)].  (2.10) 

To  simplify  this  four- dimensional  correlation  function,  we  make  an  assumption  that  the 
channel  is  wide  sense  stationary  in  the  short  term,  i.  e,  the  second-order  statistics  depends 
only  on  the  difference  between  time  instants,  rather  than  the  absolute  time.  Therefore, 
the  autocorrelation  of  h(t,  r)  is  simplified  as 

E (h(t  +  At,  ti)  •  h*(t,  r2))  =  Rh(At,  ti)<5(ti  -  r2)  (2.11) 

and  the  channel  scattering  function  is  defined  by  [11] 

/OO 

Rh(At,r)e-j2*fAtdAt  (2.12) 

-oo 

The  two-dimensional  scattering  function  (2.12)  depicts  the  power  spectral  density  of 
the  channels  along  the  geotime  and  completely  describes  the  second-order  statistics  of  a 
stationary  fading  channel. 

2.3.  Experimental  results 

The  underwater  experiment,  named  as  Reschedule  Acoustic  Communication  Ex¬ 
periment  (RACE),  was  conducted  in  Narragansett  Bay,  Rhode  Island,  by  Woods  Hole 
Oceanographic  Institution  (WHOI),  in  March  2008.  The  water  depth  varied  between  9 
to  14  meters.  The  transmitter  with  an  array  of  two  transducers  was  fixed  on  a  tripod  and 
was  four  meters  above  the  sea  bottom.  The  receiver  with  an  array  of  12  hydrophones  was 
also  fixed  on  a  tripod  and  was  located  two  meters  above  the  sea  bottom.  In  this  study, 
the  data  was  collected  when  the  distance  between  the  transmitter  and  receiver  was  1000 
meters.  The  carrier  frequency  was  fc  =  11.5  kHz,  the  sampling  rate  fs  =  39.0625  kHz, 
and  the  bandwidth  B  =  fsj 10  =  3.90625  kHz. 
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A  probing  sequence  with  a  duration  of  22  seconds  was  transmitted  to  probe  the 
channel  impulse  responses.  A  sliding  window  was  used  with  a  window  size  Np  =  128  and 
the  interval  between  adjacent  window  Nj  =  32.  Hence,  the  sampling  rate  at  the  absolute 
time  domain  is  1/(32  •  Ts)  =  122  Hz  which  means  the  channels  are  estimated  in  every  8.2 
milliseconds.  The  channel  impulse  responses  (CIRs)  were  estimated  progressively  by  the 
LS  method  described  in  (2.3).  For  each  tap,  2144  samples  were  obtained  in  one  probing 
data  packet.  The  two-dimensional  time-varying  CIRs  for  the  subchannel  corresponding 
to  (Txl,  Rxl)  and  (Tx2,  Rxl)  are  plotted  in  Fig.  2.1(a)  and  Fig.  2.1(b).  The  frequency- 
selective  channel  length  spans  25  symbol  periods. 

Fig.  2.2(a)  and  Fig.  2.2(b)  show  the  real  and  imaginary  parts  of  the  measured  data 
distribution  versus  a  Gaussian  distribution  with  the  same  mean  and  different  variances. 
We  can  see  the  measured  PDFs  of  the  real  and  imaginary  parts  do  not  follow  the  Gaussian 
distribution.  For  the  probability  densities  around  large  values,  the  Gaussian  distribution 
with  large  variance  matchs  the  tails  of  the  measured  PDFs.  However,  for  the  probability 
densities  at  small  values,  the  Gaussian  distributions  with  small  variance  fit  the  peaks  of 
the  measured  PDFs  better. 

Since  the  real  and  imaginary  parts  are  not  well  represented  by  Gaussian  distribu¬ 
tions,  it  is  reasonable  to  predict  the  magnitude  probability  distribution  will  not  follow 
a  Rayleigh  distribution.  Fig.  2.2(c)  shows  the  measured  magnitude  PDF  versus  fitted 
compound  K  and  other  distributions.  It  can  be  seen  visually  that  the  Rayleigh  and 
Gamma  distributions  are  not  a  good  fit.  The  KS  test  was  used  to  compare  the  mea¬ 
sured  PDF  with  different  distributions  to  find  the  best  match  to  the  measured  PDF. 
The  KS  test  results  for  different  distributions  are  depicted  in  Table  2.2.  In  this  table, 
the  p-value  is  associated  with  a  test  statistic,  which  represents  the  probability  of  a  test 
statistic  at  extreme  as  the  one  actually  observed.  If  the  significance  level  a  =  5%,  the 
null  hypothesis  is  rejected  when  the  p-value  is  less  than  0.05.  The  Dnn>  is  the  maximum 
distance  between  two  CDFs,  and  the  results  has  two  choice:  reject  or  accept  the  null 
hypothesis.  The  fittest  parameters  for  different  distributions  were  chosen  from  the  KS 
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10  15 

Channel  taps  (T ) 


(a)  Estimated  CIRs  for  the  subchannel  Txl- 
Rxl 


(b)  Estimated  CIRs  for  the  subchannel  Tx2 — 
Rxl 


Figure  2.1  Estimated  time- varying  frequency-selective  channel  impulse  responses. 


test  that  returned  the  highest  p- value.  The  null  hypothesis  which  asserts  the  measured 
distribution  follows  the  fitted  compound  K  distribution  with  v  =  0.63  and  /3  =  6.04  was 
accepted  with  a  p- value  of  0.96  and  Dnn>  =  0.07.  However,  the  null  hypothesis  when 
the  measured  data  is  subject  to  the  best  fit  Rayleigh  distribution  is  rejected  with  a  very 
small  p- value,  1.1  x  10~51.  The  null  hypothesis  when  compared  to  the  fitted  Gamma 
distribution  is  rejected  with  a  p- value  of  1.18  x  10-6.  Therefore,  we  conclude  that  the 
measured  data  follows  a  compound  K  distribution  rather  than  a  Rayleigh  or  Gamma 


distribution. 


PDF  PDF 
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Table  2.2 


The  Kolmogorov-Smirnov  results  for  the  measured  PDF  when  compared  to 


different  fitted  c 


istributions 


Gamma 

Rayleigh 

Compound  K 

p- value 

1.18  x  10"6 

1.1  x  KT51 

0.96 

D  / 

n,n 

0.16 

0.48 

0.07 

Results 

reject 

reject 

accept 

Normalized  real  part  Normalized  imaginary  part 


(a)  Real  Part  PDF  (b)  Imaginary  Part  PDF 

Magnitude  PDF  Phase  PDF 


(c)  Magnitude  PDF 


(d)  Phase  PDF 


Figure  2.2  The  PDFs  of  the  complex  CIR.  (a)  PDF  of  real  part  versus  Gaussian  with 
cr  =  1.6  and  0.75.  (b)  PDF  of  imaginary  part  versus  Gaussian  with  a  =  1.25 
and  0.66.  (c)  PDF  of  magnitude  versus  compound  K  and  other  distributions, 
(d)  PDF  of  phase  versus  uniform  distribution. 


The  normalized  autocorrelation  functions  of  each  channel  tap  for  real  part,  imag¬ 
inary  part  and  complex  envelop  are  depicted  in  Fig.  2.3.  It  is  noted  that  these  three 
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Figure  2.3  The  normalized  autocorrelation  of  the  real  part,  imaginary  part,  and  complex 
envelop  of  the  CIR.  The  autocorrelation  function  is  dominated  by  the  Gamma 
component  of  the  CIR. 


autocorrelation  curves  behave  similarly  all  with  exponential  decay.  Define  the  channel 
coherence  time  as  the  time  over  which  the  autocorrelation  is  above  0.5.  The  coherence 
time  was  approximately  Tc  =  0.12s  for  this  channel.  It  implies  that  the  channel  can 
be  considered  to  be  invariant  if  the  length  of  data  block  is  less  than  Tc  ■  B  ps  470  sym¬ 
bols.  The  relation  between  the  coherence  time  and  the  maximum  Doppler  spread  is  given 


by  [12] 


Tc 


9 

167T  fm 


(2.13) 


where  fm  is  maximum  Doppler  shift.  In  this  case,  fm  =  1.5  Hz  and  the  dynamic  range  of 
the  instantaneous  Doppler  spread  is  [—1.5, 1.5]  Hz.  Figure  2.4  shows  the  cross  correlation 
of  the  real  and  imaginary  part  of  the  CIR.  It  is  seen  that  the  real  parts  and  the  imaginary 
part  are  uncorrelated. 

The  two-dimensional  scattering  function  is  shown  in  Fig.  2.5  by  performing  Fourier 
transform  on  the  correlation  function.  We  see  that  the  average  Doppler  shift  is  very  small, 
close  to  zero  Hz  since  the  transceiver  platforms  were  always  fixed  during  the  experiment 
and  the  motion  of  the  wave  was  not  significant.  From  the  power  spectral  density,  the 
channel  has  a  Doppler  spread  of  [—1.5, 1.5]  Hz  at  a  delay  of  1  ~  2.5  ms. 


30 


Crosscorrelation  of  real  and  imaginary  parts 


Figure  2.4  The  normalized  cross  correlation  between  the  real  and  imaginary  parts  of  the 
CIR.  The  real  and  imaginary  parts  are  uncorrelated. 


Scattering  function 
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Doppler  Spread  (Hz) 


Figure  2.5  The  channel  scattering  function 


2.4.  Conclusion 

In  this  paper,  we  presented  the  statistical  characteristics,  including  PDF  and  second- 
order  statistics,  for  the  UWA  channels  measured  in  RACE,  2008.  Experimental  results 
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demonstrate  that  the  real  and  imaginary  parts  of  the  complex  UWA  channels  do  not  fol¬ 
low  Gaussian;  the  PDF  of  the  magnitude  fits  the  compound  K  distribution;  the  PDF  of 
phase  approximately  follows  uniform  distribution.  The  autocorrelations  of  CIRs  behave 
like  decaying  exponential  curves  and  the  channel  coherence  time  is  also  roughly  esti¬ 
mated.  The  scattering  function  provides  better  understanding  on  the  spectral  density  of 
the  UWA  channels. 
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Abstract 

The  spatial  and  intertap  correlation  matrices  of  fading  channels  is  first  estimated 
for  underwater  multiple  input  multiple  output  (MIMO)  systems  using  real  experimental 
data.  We  confirm  that  underwater  MIMO  channels  exhibit  strong  spatial  and  temporal 
correlation.  The  bit  error  rate  (BER)  performance  of  the  MIMO  system  at  different 
correlation  scenarios  are  also  evaluated  with  computer  simulation.  At  BER=0.01%,  the 
required  SNR  of  highly  correlated  channel  is  5  dB  higher  than  the  i.i.d.  channel.  BER 
results  of  the  experimental  data  also  quantitatively  show  the  performance  degradation 
caused  by  the  correlation. 

3.1.  Introduction 

Underwater  acoustic  (UWA)  transmission  is  an  effective  method  used  for  commu¬ 
nication  over  medium  and  long  ranges  (1-1000  km);  however,  many  aspects  of  UWA 
channels  present  additional  challenges  when  compared  to  in  air  radio  frequency  (RF) 
channels.  UWA  channels  exhibit  time  variance,  large  multipath  delay  times,  frequency 
dependent  fading,  and  slow  wave  propagation  speed.  The  time  variance  is  caused  by 
changes  in  the  environment  such  as  changes  in  the  speed  of  the  transmitters,  receivers, 
or  the  media, itself,  and  turbulence  caused  by  surface  waves  and  vehicles.  Additional  com¬ 
putational  complexity  results  when  the  channel  must  be  re-estimated,  repeatedly.  The 
multipath  delay  times  may  be  on  the  order  of  tens  to  hundreds  of  milliseconds  resulting 
in  great  intersymbol  interference  (ISI)  and  channel  impulse  responses  (CIR)  that  have 
many  more  taps  than  an  RF  CIR.  The  frequency  dependent  fading  severely  limits  the 
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bandwidth  available  for  data  transmission.  Waves  that  propagate  more  slowly  are  more 
affected  by  small  Doppler  shifts  and  spread. 

Data  collected  in  the  Surface  Process  and  Acoustic  Communications  Experiment 
conducted  in  fall  2008  (SPACE08)  by  the  Woods  Hole  Oceanographic  Institute  (WHOI) 
is  investigated.  The  effects  of  spatial  correlation  on  the  BER  are  investigated.  It  is  shown 
that  the  a  smaller  amount  of  receivers  may  be  used  in  place  of  a  larger  number  of  receivers 
with  a  very  small  impact  on  the  BER  as  long  as  the  aperture  is  the  same.  The  Kronecker 
properties  of  the  spatial  correlation  matrix  are  analyzed.  The  matrix  is  decomposed  into 
a  series  of  Kronecker  products  and  the  error  resultant  matrix  formed  by  the  decomposed 
matrices  and  the  original  matrices  is  quantified.  A  simulation  is  conducted  to  confirm 
the  experimental  results  are  due  to  spatial  correlation. 

3.2.  MIMO  Acoustic  Communication  Systems 

The  system  may  be  generalized  as  a  MIMO  system  with  P  transmitters  and  Q 
receivers.  The  baseband  equivalent  discrete  time  domain  signal  at  the  qth  receiver  is 
defined  as 

p  L 

yq(k)  =  5Z  5Z  k)xp(k  +  1  -  0  +  wq(k)  (3.1) 

p= l  i=i 

where  L  is  the  channel  length;  hq:P(l,  k)  is  the  lth  CIR  coefficient  at  time,  k ;  T  is  the 
symbol  interval;  fq,p,k  is  the  instantaneous  Doppler  shift  between  transmitter,  p,  and 
receiver,  q ;  6qp  is  the  initial  phase  error  present  in  the  q,p  subchannel;  and  wq(k)  is  zero 
mean  white  Gaussian  noise  with  variance,  present  at  receiver,  q. 

The  complete  MIMO  system  is  modeled  in  matrix  form  as 


hi,i  •  • 

^l.P 

Xl 

Wi 

y  = 

^Q.l 

k Q,P 

.  Xp  . 

+ 

w 

Q 

=  H  x  + w 


(3.2) 
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where  the  received  vector  is  y  =  [yf,  •  •  •  ,  y7)]7  .  The  signal  received  at  each  hydrophone 
may  be  represented  as  yq  =  [yq(  1),  •  •  •  ,  yq(N)]T.  xp  =  [xp(l),  •  •  •  ,xp(N)}T  is  the  signal 
transmitted  by  the  pth  transducer  where  xp([Nd  +  1,  •  •  •  ,  TV])  =  0. 

=  dia,g{ei(2nTfq’p’1+eq’p\  ■  ■  ■  ,  e^2nNT^'p’N+dq’p^ .  hQyP  is  the  ( p,q)th  CIR  with  the 
diagonal  phase  components  missing. 

It  is  assumed  the  coherence  time  of  the  channel  is  greater  than  the  block  duration 
which  implies  the  channel  may  be  considered  time  invariant  over  one  block. 

3.2.1.  Channel  Estimation. 

A  time  domain  least  squares  method  is  used  to  estimate  the  channel  from  the  data 
contained  in  the  transmitted  blocks.  The  transmitted  data  matrix  at  the  pth  transmitter, 
xtd,  is  defined  as 


Xtd(L) 

xtd{L  -  1)  •  ■ 

Xtd(  1) 

xtd(L  +  1) 

xtd(L) 

Xtd{  2) 

Xtd(Nd ) 

xtd(Nd  -  1)  •  • 

xtd(Nd— L+l) 

(3.3) 


where  Nd  is  the  number  of  data  symbols  used  in  the  estimation. 

The  subchannel  impulses  involving  receiver,  q ,  are  defined  as  hg  =  [A^ih7  1;  •  •  •  ,  \q^ph^p]T . 
The  estimation  is  performed  using  (3.4). 


=  xj  •  ytq  (3.4) 

where  xt  =  [xn,  •  •  •  ,  xt] ,  (•)t  denotes  the  pseudo- inverse,  and  ytq  =  [ytq(L),  •  •  •  ,  ytq(Np)]T 
denotes  the  received  data  symbols  at  hydrophone,  q. 

3.2.2.  Frequency  Domain  Turbo  Equalization. 

The  incoming  data  blocks  are  converted  to  the  frequency  domain  with  a  simple 
FFT.  The  vector  Yp  is  the  frequency  domain  representation  on  the  incoming  data.  The 
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MMSE  FDTE  adaptively  computes  new  coefficients  using  the  a  priori  mean  and  variance 
of  the  incoming  data  in  order  to  estimate  the  symbols  sent. 

8PSK  was  used  as  the  modulation  technique  of  all  estimations.  The  symbol  alphabet 
used  is  5  =  [1,  ^=(1  +j),j,  ^(-1  +j),  -1,  ~j,  7^(1 -/)]  for  the  bit  patterns, 

Ki,i4,24i,3]  €  {111,  110,  010,  000,  lOO,  101,  OOl,  Oil},  respectively.  The  mean,  /ipk,  and 
variance,  up^,  are  calculate  using 

=  1/4  •  ((/ 1  —  1)  •  (I3  +  I2 )  —  V2  ■  {I2  +  h  ■  I3 )) 

+j  •  1/4  •  ((/1  +  1)  •  (Z3  —  I2)  +  V2  ■  (I3  +  /1  •  I2 )) 

1  |hp,fc|  (3.5) 

where  lg  =  tanh(l//(cP)3(fc_i)+s)/2).  cpy  is  the  coded  bit  and  !//(•)  is  the  a  prion  log 
likelihood  ratio  of  the  decoded  bit. 

The  MMSE  equalized  symbols  are  represented  in  the  frequency  domain  as 


X,(fcj  =  Kp  •  U"  •  (Y  -  H  •  X  +  ,.,,IIF,uM) 


(3.6) 


where  H  is  the  frequency  domain  channel  matrix;  F p  is  the  normalized  DFT  matrix;  and 
uPifc  is  a  vector  of  length  PN  where  element,  ((p—l)N  +  k),  is  1  and  the  other  elements 
are  zero.  The  vector  X  is  the  mean  of  the  incoming  symbols  in  frequency  domain  as  is 
defined  as 


X 


*■1  j 


>xpJ 


T-iT 


(3.7) 


where  xp  =  \jiPt i,  •  •  -  ,  hp,n]t ■ 

The  equalizer  coefficients,  Up,  are 


H , 


-1 


H 


O 


O 


( p-l)NxN 

In 

( P-p)NxN 


(3,8) 
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where  Omxn  is  a  zero  matrix  and  V  is 


V  —  diag{i'1IN, 


(3.9) 


where  vp  is  the  mean  of  up  and  the  scalar  Kp  is 

1  V 

Kp  =  1  H  Tr  {H  at_  tv+  j  :pAr)  Up } ,  (3.10) 

where  denotes  the  matrix  composed  of  the  i-th  to  j-th  columns  of  H. 

The  time  domain  symbols  for  the  transmitter  in  (3.6)  are 

xp(k)  =  k)  ■  Xp(fc),  (3.11) 


where  F(#  k)  represents  column,  k  of  the  DFT  matrix  F. 

3.2.3.  Spatial  Correlation. 

A  MIMO  system  may  be  modeled  as  Y  =  HX  +  W  where  channel  matrix,  H,  is 
defined  by 


H  = 


(3.12) 


h 


04 


h 


Q,p 


In  [1],  the  spatial  correlation  is  defined  as  the  Kronecker  product  of  the  transmitter, 
receiver,  and  intertap  correlation  matrices.  The  spatial  correlation  between  receivers  q\ 
and  q-2  is  defined  as 


Rqiq2  =  E{hvechTvec}  =  <g>  ®  Ttap  (3.13) 

where  in  a  two  transmitter  MIMO  system,  hvec  =  [hqipi,  hqiP2,  hqnn ,  hq2P2\.  The  Kronecker 
product  is  represented  by  <g).  Matrices  Ttx,  'hfe,  and  Ttaf>,  are  the  transmitter,  receiver, 
and  intertap  correlation  matrices,  respectively. 
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The  transmitter  correlation  matrix  may  be  factored  out  by  using 

Rgig2  =  E{hRxhTRx}  =  T Rx  <g>  T tap 

^dld^tap  ^qiq2^tap 
^9291  ^  tap  ^9292  ^  tap 

where  hRx  =  [hqipi,  hq2P1). 

A  metric  that  may  be  used  to  quantify  the  similarity  of  two  correlation  matrices 
was  described  in  [2],  The  correlation  matrix  distance  (CMD)  is  defined  as 

dcorr(  Rl!R2)  =  l-^2^4-  (3-15) 

ii^mI  i/I  ik2||/ 

where  tr{-}  is  the  trace  operator  and  ||  •  ||/  is  the  Frobenius  norm.  The  CMD  will  return 
as  zero  for  matrices  that  are  identical  save  for  a  scaling  factor  and  one  for  matrices  that 
are  maximally  distant. 

The  CMD  is  based  upon  the  inner  product  and  measures  the  angle  between  the 
two  correlation  matrices  in  n  x  n  dimensional  space.  Matrices  that  are  identical  up  to  a 
scaling  factor  have  no  angle  between  them  so  the  CMD  will  return  1  —  cos(O)  =  0. 

The  CMD  is  useful  when  comparing  a  measured  correlation  matrix  to  an  estimated 
correlation  matrix  because  it  gives  a  quantifiable  method  to  determine  the  difference 
between  the  spatial  structure  (and  therefore  channel  characteristics)  of  the  measured 
and  estimated  channels. 

3.3.  Experiment  Results 

This  paper  analyzes  channel  spatial  correlation  through  data  collected  in  the  Surface 
Process  and  Acoustic  Communications  Experiment  conducted  in  fall  2008  (SPACE08) 
by  the  Woods  Hole  Oceanographic  Institute  (WHOI).  The  experiment  took  place  in 
shallow  water  (about  15  meters  deep)  south  of  Cape  Cod,  MA.  The  distances  between 
the  transmitter  array  and  receiver  array  were  200  m  and  1000  m.  The  transmitters  were 


(3.14) 
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Figure  3.1  Magnitude  of  the  channel  correlation  matrix  of  a  2x2  MIMO  channel 


fixed  on  a  tripod  with  the  top  transmitter  3  m  above  the  sea  floor  with  50  cm  spacing 
between  the  transmitters.  The  receivers  were  also  fixed  upon  a  tripod.  At  200  m,  there 
were  24  receivers  with  a  spacing  of  5  cm.  At  1000  m,  there  were  12  receivers  with  a 
spacing  of  12  cm.  The  carrier  frequency  was  13  kHz;  the  sampling  frequency,  Fs,  was 
39.0625  kHz;  and  the  bandwidth  was  Fj,  =  Fs/A  =  9.7656  kHz. 

3.3.1.  Spatial  Correlation. 

Measurements  taken  during  the  experiment  show  a  decrease  in  spatial  correlation 
as  the  receivers  are  spaced  further  apart.  Figure  3.1  shows  the  spatial  correlation  matrix 
between  receivers  one  and  two.  The  correlation  matrices  between  other  receivers  look 
similar.  The  shape  of  the  matrix  strongly  suggests  the  Kronecker  based  model  in  [1]  is 
applicable. 

The  decrease  in  spatial  correlation  may  be  observed  when  using  (3.14)  to  factor  out 
the  transmitter  correlation.  Figure  3.2  shows  the  spatial  correlation  between  receivers 
1  and  2,  1  and  4,  and  1  and  10.  This  figure  shows  the  maximum  and  minimum  spatial 
correlation  present  in  the  different  receiver  combinations  analyzed  above,  ft  is  evident 
that  the  spatial  correlation  decreases  with  distance. 
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(a)  ^2,1  ^tap  (b)  ^4,1  ^tap 


x  10*3 


(C)  4'10,1’I'tap 

Figure  3.2  Spatial  correlations  between  different  receivers 


Equation  (3.14)  may  be  used  to  factor  4/tx  out  of  (3.13).  By  using  the  definition  of 
the  Kronecker  product  in  reverse  and  normalizing  the  correlations  matrix  in  respect  to 
the  highest  peak  in  the  overall  correlation  matrix,  the  transmitter  and  receiver  correlation 
matrices  may  be  may  be  factored  out.  Representative  transmitter  and  receiver  matrices 
for  the  spatial  correlation  between  different  receivers  are  shown  in  Fig.  3.1. 

The  CMD  calculated  using  the  estimated  matrices  in  Fable  3.1  and  (3.15)  are  0.24, 
.17,  and  .25,  respectively.  It  is  evident  from  these  correlation  matrices  that  the  spatial 
correlation  is  decreasing  as  the  receivers  move  further  apart. 
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Table  3.1  d+x  and  T^x  f°r  hydrophones  at  200  m 


and 

h2 

dbx  = 

dRx  = 

'  1.0000  +  jO.OOOO  0.8861  -  j0.1140  ' 
0.8861  +  j'0.1140  0.7983  + jO.OOOO 

'  0.7606  +  jO.OOOO  0.7692  +  j0.4101  ' 
0.7692  - j0.4947  1.0000  +  jO.OOOO 

and 

H± 

dbx  = 

d'rx  = 

'  1.0000  +  jO.OOOO  0.8954  -  j0.2687  ' 
0.8954  + j0.2687  0.8741  +  jO.OOOO 

"  0.6090  +  jO.OOOO  0.7158  +  j0.3098  ' 
0.7158- j0.3098  1.0000  +  jO.OOOO 

H, 

and 

d'tx  = 

'  0.8551  +  jO.OOOO  0.7077  -  j0.5951  ' 
0.7077  + j0.5951  1.0000  +  jO.OOOO 

Hw 

d'rx  = 

'  0.4018  -  jO.OOOO  0.5961  +  j0.2142  ' 
0.5961  -j0.2142  1.0000  +  jO.OOOO 

3.3.2.  Transceiver  BER  Performance. 

The  BER  of  MIMO  configurations  with  different  receivers  spaces  was  calculated.  It 
was  found  that  a  larger  space  between  the  receivers  had  smaller  spatial  correlation  and 
therefore  lower  BERs.  The  BERs  of  60  runs  of  200  m  data  and  24  runs  of  1000  m  data 
conducted  during  SPACE08  were  calculated  while  using  a  2  iteration  turbo  equalization 
technique.  All  runs  analyzed  used  8PSK  modulation.  First,  the  average  BER  was  cal¬ 
culated  for  receivers  1,  2,  3,  and  4.  Then  the  average  BER  of  receivers  1,  4,  7,  and  10, 
was  calculated.  Finally  the  average  BER  using  all  receivers  from  1  to  10  was  calculated. 
The  average  BERs  were  calculated  for  200  m  and  1000  m  data.  The  results  are  shown 
in  Table  3.2. 

The  MIMO  receiver  system  with  receivers  1,4,7,  and  10,  has  the  same  aperture  as 
a  system  with  receivers  1  through  10  but  less  receiver  power  power  due  to  the  smaller 
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Table  3.2  Transmission  distance  and  receiver  spacing  BER  results 


1,2, 3, 4 

1,4,7,10 

1  through  10 

200  m 

1.58e-l 

3.03e-2 

9.97e-3 

1000  m 

2.23e-l 

9.3E-2 

9.29e-2 

amount  of  receivers.  The  MIMO  system  using  receivers  1,2,3,  and  4,  has  the  same 
receiver  power  as  the  system  using  receivers  1,4,7,  and  10,  yet  it  has  a  smaller  aperture. 
The  results  show  that  the  aperture  is  more  important  than  receiver  power.  The  slightly 
smaller  BER  of  the  ten  transmitter  system  is  explained  by  the  greater  receiver  power; 
however,  the  ten  transmitter  system  is  not  as  efficient  as  the  system  with  1,4,7,  and 
10,  receiver  spacing.  At  200  m,  it  took  six  extra  receivers  to  decrease  the  1,4,7,10  BER 
by  about  37  percent,  yet  merely  spacing  out  the  same  number  of  receivers  was  able  to 
decrease  the  1,2, 3, 4  BER  by  about  97  percent.  At  1000  m,  the  advantage  of  increased 
distance  between  receivers  is  more  evident.  The  BERs  between  all  ten  receivers  and  the 
4  spaced  receivers  are  nearly  identical. 

3.4.  Simulation  Results 

We  conducted  computer  simulation  to  verify  the  effect  of  channel  spatial  correla¬ 
tion  on  BER  performance  of  the  FDTE  receiver.  In  the  simulation,  a  MIMO  2x4  UWA 
system  was  considered  with  80-tap  frequency-selective  channels.  The  maximum  Doppler 
spread  was  fd  =  1  Hz  and  the  symbol  period  was  T  =  0.1024  ms.  We  adopted  two  sparse 
channel  Power  Delay  Profiles  (PDP),  as  shown  in  Fig.  3.3(a),  where  the  total  average 
power  of  each  profile  was  normalized  to  one.  The  spatial-temporally  correlated  fading 
channel  waveforms  were  generated  by  the  triply-selective  simulation  model  [1]  with  each 
ray  following  Rayleigh  distribution  [7].  The  overall  channel  distribution  of  80  taps  was 
close  to  the  compound  K  distribution  as  indicated  by  recent  research  on  UWA  fading 
channels  [4,5].  The  spatial  correlation  was  simulated  via  transmit  and  receive  spatial  cor¬ 
relation  matrices,  whose  elements  were  denoted  4 ’tx(Pi,Pj)  and  vhr:r(<7,  Qj),  respectively. 
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(a)  Sparse  Channel  PDP 


n - 1 - T - 1 - T - 

— O  PDP1 
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Figure  3.3  Simulation  results:  (a)  PDP  of  the  UWA  channels  with  80  taps.  The  normal¬ 
ized  maximum  Doppler  spread  was  /^Ts  =  1.024  x  10”4.  (b)  BER  performance 
of  the  FDTE  receiver  on  different  spatial  correlations. 


We  selected  'Prx(w  (-lj)  —  to  approximate  the  spatial  correlation  among  the  hy¬ 

drophone  elements  since  the  hydrophones  were  equally  spaced  in  the  experiment.  The 
encoded  bits  were  interleaved  randomly  and  were  mapped  into  8PSK  symbols  via  Gray 
coding. 

We  applied  the  FDTE  algorithm  to  recover  the  transmit  data  streams  and  computed 
the  average  BER  via  10  to  70  trials  with  76,800  information  bits  per  trial.  The  BER 
performance  in  channels  of  different  spatial  correlations  is  depicted  in  Fig.  3.3(b).  The 
BER  performance  in  the  i.i.d  (independent,  identically  distributed)  fading  channels  is 
also  included  in  Fig.  3.3(b)  as  bench  marks. 

The  spatial  or  intertap  correlation  of  the  MIMO  channels  had  strong  impact  on 
the  receiver  performance,  as  illustrated  in  Fig.  3. 3(b). The  i.i.d.  channel  had  no  spatial 
nor  intertap  correlation;  while  the  channel  with  xVtx  —  0  and  T rx  =  0  had  no  spatial 
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correlation  but  intertap  correlation.  The  performance  gap  between  these  two  channels 
was  1.1  dB  at  BER=10-4.  As  the  transmit  correlation  or  receive  correlation  increases 
on  top  of  the  intertap  correlation,  the  BER  performance  becomes  worse.  At  the  extreme 
case  of  Tfe,,  =  0.9  and  T rx  =  0.9,  the  BER  performance  was  about  4  dB  worse  than  that 
of  no  spatial  correlation,  or  5  dB  worse  than  that  of  the  i.i.d.  channel. 

3.5.  Conclusion 

In  this  paper  the  effects  of  spatial  correlation  on  the  BER  are  analyzed.  The  results 
show  that  the  BER  is  relatively  unaffected  by  a  reduction  in  receivers  as  long  as  the 
aperture  is  the  same.  A  correlation  matrix  was  decomposed  into  an  estimated  Kronecker 
product  and  the  error  between  the  original  matrix  and  the  estimated  Kronecker  product 
was  quantified.  Simulations  with  differing  amounts  of  spatial  and  iter-tap  correlation 
were  conducted.  The  simulation  and  experimental  results  show  that  the  BER  of  MIMO 
systems  with  high  spatial  correlation  is  worse  than  systems  with  low  spatial  correlation. 
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SECTION 


2.  CONCLUSIONS 


Statistical  characteristics,  including  PDF  and  second-order  statistics,  for  the  UWA 
channels  measured  from  experimental  data  are  analyzed.  Experimental  results  demon¬ 
strate  that  the  real  and  imaginary  parts  of  the  complex  UWA  channels  do  not  follow 
Gaussian;  the  PDF  of  the  magnitude  fits  the  compound  K  distribution;  the  PDF  of 
phase  approximately  follows  uniform  distribution.  The  autocorrelations  of  CIRs  behave 
like  decaying  exponential  curves  and  the  channel  coherence  time  is  also  roughly  esti¬ 
mated.  The  scattering  function  provides  better  understanding  on  the  spectral  density  of 
the  UWA  channels.  Effects  of  spatial  correlation  on  the  BER  are  also  analyzed.  It  is 
shown  that  the  BER  is  relatively  unaffected  by  a  reduction  in  receivers  as  long  as  the 
aperture  is  the  same.  Experimental  results  show  that  the  BER  of  MIMO  systems  with 
high  spatial  correlation  is  worse  than  systems  with  low  spatial  correlation.  Simulations 
with  differing  amounts  of  spatial  and  iter-tap  correlation  were  conducted  verifying  the 
experimental  results. 
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